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Abstract 

The formalism which has been developed to give general expressions for the 
determinants of differential operators is extended to the physically interesting 
situation where these operators have a zero mode which has been extracted. 
In the approach adopted here, this mode is removed by a novel regularisation 
procedure, which allows remarkably simple expressions for these determinants 
to be derived. 
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I. INTRODUCTION 



The increasing use of path integrals as a calculational tool has led to a corresponding 
increase in interest in the evaluation of functional determinants. This is simply because the 
evaluation of Gaussian path integrals typically gives such determinants. The first results 
were obtained over thirty years ago: Gel'fand and Yaglom JIJ derived expressions for the 
functional determinants obtained from evaluating path integrals with the simplest type of 
quadratic action. In subsequent years the results became more general and the formalism 
more elaborate culminating with the work of Forman || who gave a remarkably simple 
prescription which can be applied to a rather general operator and boundary conditions. 

However, in many calculations involving Gaussian integrals which are currently carried 
out, these results are not directly applicable. The reason is that the Gaussian nature of the 
integral is frequently a consequence of expanding about some non-trivial "classical" solution 
of the model (e.g. a soliton or instanton). Typically this results in a particular point (in 
space or time) being selected, which breaks the translational invariance of the theory, and so 
gives rise to a Goldstone mode. There are other possible ways that such a zero mode could 
come about, but in all cases the Gaussian approximation breaks down. The remedy is to first 
extract this mode as a collective coordinate || and to treat only the non-zero modes in the 
Gaussian approximation. Therefore, it is not the functional determinant which is required 
in these cases — it will in any case be identically zero — but the functional determinant 
with the zero mode extracted. 

In this paper we present a systematic method to calculate this quantity. The most 
obvious way to proceed is to "regularise" the theory in some way, so that the eigenvalue of the 
operator under consideration which was previously zero is now non-zero. The determinant 
is now also non-zero and the pseudo-zero eigenvalue can be factored out, the regularisation 
removed, and a finite result obtained. Previous approaches have been rather ad hoc, being 
performed on a case by case basis as the need arose. For example, it may be possible 
in certain cases to modify the form of the operator in such a way that the zero mode is 
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regularised, but also that the calculation may still be performed 0,0). In other cases, it 
may be possible to move the boundaries to achieve the same end |8],|9|]. Here we adopt an 
approach which applies to very general situations and which, we believe, is the simplest and 
most systematic regularisation and calculational procedure. This is because the method is 
the least intrusive — the operator and the position of the boundaries are left unchanged 
- only the form of the boundary conditions are modified in the regularisation procedure. 
We use the notation and general approach of Forman to calculate the regularised functional 
determinant, since it is ideally suited to this form of regularisation, emphasising as it does 
the separation of the boundary conditions from the solutions of a homogeneous differential 
equation. 

The outline of the paper is as follows. In section 2 we develop our method in one of 
the simplest situations in order to clearly illustrate it. The calculation of the regularised 
expression for the formerly zero eigenvalue is derived, for the most general case that will 
interest us, in section 3 and the general procedure for finding the functional determinant 
with the zero mode extracted is described in section 4. In section 5 we apply the method to 
certain specific cases and we conclude in section 6 with some general remarks. 



II. A SIMPLE EXAMPLE 

In this section we will explain the method by carrying out an explicit calculation on what 
is perhaps the simplest example. Suppose that we wish to calculate the determinant of an 
operator of the form 

where P(t) is a known real function. We suppose that the boundary conditions on the 
functions on which L operates is u(a) = u(b) = 0. In particular, the eigenfunctions of L 
have to satisfy these conditions. 

We now give Forman's prescription for calculating det L. A more detailed discussion is 
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given in section 4, where our approach is explained in greater generality. The recipe has two 
ingredients: 

(i) Write the boundary conditions on L in the form 

"01 

(2.2) 

where M and N are 2x2 matrices and u = du/dt. These two matrices are not unique; for 
the case of our boundary conditions u(a) = u(b) = we choose them to be 
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(2.4) 



(2.3) 

(ii) Now consider a different problem. Let yi{t) and y 2 {t) be two independent solutions 
of the homogeneous differential equation Lh = 0. Construct 

'yi{t) Mt)' 

and the 2 x 2 matrix Y(b) = H(b)H- l (a). 

Forman then proves that Q 

detL _ det(M + NY(b)) 
detL ~ det(M + A^F(6)) 

We would expect that detL itself is divergent, being a product of an infinite number of 

eigenvalues of increasing magnitude. Therefore it is only when it is defined relative to the 

determinant of an operator of a similar type (denoted here by L), that it has any meaning. 

In applications to path integrals, ratios of determinants such as the one on the left-hand-side 

of ( |2.5| ) naturally arise from the normalisation of the path integral itself. In general, they 

will relate to a simple quantum system or stochastic process, such as the harmonic oscillator 

or Ornstein-Uhlenbeck process. In these cases, P(t) is independent of t and will not, in 

general, have a zero mode. 

For the matrices M and TV of our simple example, 



det(M + iVy(&)) = Y 12 (b) 

= yi(a)y2(b) -y 2 (a)yi(b) 
yi(a)y 2 (a) - y 2 (a)y 1 (a) 



(2.6) 



The denominator of this expression is the Wronskian, which does not vanish since the two 
solutions yx{t) and y 2 (t) are presumed independent. If we take y\{t) to be a solution for 
which yi(a) = 0, then ( p76| ) can be simplified to yi(b)/yi(a), so that, if yi(a) = too, 

detL = z/i(6)yi(6) ? , 
detL Vi(b)yi(b) 

This simple expression is particularly useful, since it only involves y\,y\ and their first 
derivatives at one of the boundaries. We should stress that results such as this have been 
known since the work of Gel'fand and Yaglom — our purpose here is to introduce the 
formalism required to describe our approach, in as simple a way as possible. 

Now suppose that y\{b) = (as well as y\{a) = 0). Then yi(t) is an eigenvalue of L with 
zero eigenvalue. This is the situation of interest to us in this paper. To extract this zero 
mode, we first regularise the problem by modifying it so that the operator is unchanged, 
but the boundary conditions u(a) = u(b) = become 



w (e) (a) = 0, u {e) (b) = eu ( '\b) 



(2.8) 



where e is some small number. So now yi(t) is no longer an eigenfunction of L with zero 
eigenvalue. Let y[ e \t) be the corresponding eigenfunction (i.e. the one which reduces to yi(t) 
when e — ^ 0) and let it have eigenvalue \^ e \ To find detL with these boundary conditions 
we first note that Y(b) is unchanged, since it does not involve boundary conditions at all; it 
only depends on two independent solutions to the homogeneous differential equation Lh = 0. 
Modifying the boundary conditions as in (|2.8| ) only changes M and iV to 



"1 0" 




"0 " 








0_ 




1 -e_ 



This gives det(M( e ) + N^Y(b)) = Y 12 (b) - eY 22 (b). But since Vl {a) 
and so 



(2.9) 

yi(b) = 0, Y 12 (b) = 0, 



det(M (e) +N (e) Y(b)) 



-tY 22 {b) 
Mb) 



(2.10) 



This is the regularised form of the determinant. In the next section we will give a general 
method for finding A^. In this simple problem it turns out that, to lowest order, 

A (e) = _ e /lW (2.11) 

{vim) 

where (yi\yi) is the norm of the zero mode: 



{Vi\Vi)= / dty((t). (2.12) 

From (|2.10| ) and (|2.11|) we have that 

lim det(MM+OT)) = (yM 
e-o yx{a)yi{b) 

This remarkably simple result is the one which we sought. Note that, apart from the norm, 
it is only involves y\ at the boundaries. In applications, it will usually be the case that 
the norm in ( |2.13p will cancel with an identical factor coming from the lowest order form 
of the Jacobian of the transformation to collective coordinates. Therefore (y\\yi) need not 
be calculated. Denoting the determinant of L with the zero mode extracted by det' L and 
normalising by det L, we finally obtain 

det'L = ( yi \ yi ) y^b) 
detL yi(a)yi(b) yx(b) 

The method we have described to find the regularised form of det(M + NY (b)) is hardly 
more complicated than finding the unregularised form. The key to achieving this happy state 
of affairs was firstly the decision to modify only the boundary conditions, and secondly, the 
choice of regularised boundary conditions which gave simple forms for M^ € > and We 
shall now show that these choices also allow \^ to be determined in a very simple and 
elegant way. 



III. THE REGULARISATION OF THE EIGENVALUE 



While the regularised form of the determinant could be found by use of Forman's method, 
a new technique for calculating the previously vanishing eigenvalue, \( e \ has to be developed. 
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It is natural to attempt to calculate it perturbatively in e, but it is not at all obvious that a 
general procedure can be set up. Fortunately, it will turn out that choosing the regularised 
boundary conditions in the manner illustrated in section 2 on a simple example, enables 
to be found to lowest order almost without calculation. 

Let us begin describing the method where the operator is of the simple form ( |2.1p ; we will 
generalise to more complicated operators later in this section. There is no need to specify 
the boundary conditions at this stage, since, as we will see, a useful formula for can 
be derived without having to make any choices of boundary conditions. Using the notation 
introduced in the last section 

b$ = \^y[ e) (3-1) 
where y[ e \t) — > yi(t) and A^ — > as e — > 0. From ([3.1|) 

[ b dt yi Lyi e) = \U f b dt yi y[ e) 

J a J a 

= ^ ] (yi\yi) , (3.2) 



to lowest order in e. Integrating by parts gives, again to leading order, 



yi e \t)y 1 (t)-y 1 (t)yi e) (t)] b 



(3-3) 



(Vi\Vi) 

This result is true for operators of the form (|2.1| ) with arbitrary boundary conditions. As an 
example, suppose we impose the regularised boundary conditions ( |2.8| ). Then the eigenfunc- 
tion y[ e \t) will satisfy them: y[ (a) = 0, y[ e \b) = ey[ e \b). In addition y±(a) = yi(b) = 0, so 
that to lowest order 



{yM 
yi{b)yi{b) 



(3.4) 



<2/i|z/i> 

as given in section 2. Note that the e dependence simply came from the requirement that 
y{ e \b) = eyi(b), to lowest order. 

Analogous results to ( jOD hold for more general operators. For example, suppose that 



d 2 



(3.5) 



where P(t) is a complex matrix, and suppose that the operator (|3.5|) has a single zero 
mode yi t i, that is, X^=i ^ijVjA = 0- I* 1 matrix notation, the zero mode is the column 
vector y = (2/1,1, . . . , y r ,i) T ■ Let y^(t) be the corresponding eigenfunction of the regularised 
problem with eigenvalue A^. Then to lowest order 

rb 



f dtJ2ylL ijy ^ = X {e) j:(yi tl \y^) 

J a 



(3.6) 



where now 



(l/ill/i) = E(^il^i) = ["(It £ |^,i(t)| 2 (3.7) 
Integrating the left-hand-side of ( |3.6|) by parts gives the leading order result 

ELi - ymv^t , X> v,,//;;,(/){/^ - Pfty$(t) 



+ 



(3.8) 



<2/i|l/i) (2/1I2/1) 
In most cases of interest to us L will be formally self-adjoint, and so the second term in ( p.8|) 
will vanish. The self-adjoint nature of L is expected from its origin as the second functional 
derivative of the action in the path integral with respect to the fields: 
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L(t,t') 



(3.9) 



H3 bu\(t)bUj(t') 

The most general operator which we will study in this paper takes the form 

Lij = [fl>(*)]«^ + i P ^)hj t + [P*(t)]v (3.10) 



where Po(t), Pi(t) and P2{t) are complex r x r matrices. We begin by making the transfor- 
mation 

p ii (t)=exp{i/ t t Zt(P )- 1 (Pi)} , Pi j (t)=\p{P )-\P 2 )(p)- 1 l.-\p(p)- 1 ].. (3.11) 
so that 



Lij = {P )ik ijp X )kl£>lm ip) 



rnj 



(3.12) 



where 



lJ dt 2 

Now if L is such that Po(t) = Po(t), then 

det L det C 



Ai = ^4+^«(*) (3-13) 



(3.14) 



det L det £ 

where £ is as in ( [3.13|) , but with P replaced by P. Therefore the problem has been reduced 
to that considered earlier in this section (see ( |3.5|) et seq). In fact, as regards determining the 
ratio of the determinants, Forman gives a general expression for the left-hand-side of (|3.14 ) 
(see next section), and so there is no need to implement the transformation ( ft.lip . To find 
the eigenvalue \( € \ however, this transformation is useful. It is easy to see that C has a zero 
mode if, and only if, L does, and that, in particular, if yi{t) is an eigenfunction of with 
zero eigenvalue, then Zi(t) = J2jPij{t)yj(t) is an eigenfunction of with zero eigenvalue. 
The results ( |3.6| - |3.8| ) now hold, but with L and y replaced by L and z, respectively. As in 
all of the examples discussed in this section, a judicious choice for the boundary conditions 
on the regularised eigenfunction y^' will yield an explicit regularised form for A^ with the 
minimum of calculational effort. 



IV. GENERAL PROCEDURE 

There are two aspects to our approach to the calculation of det'L/ detL. One is the 
operation of finding A^ to leading order, which was explored for the general case in the last 
section. The other aspect concerns the application of Forman's method for the calculation of 
det L/ det L, but with the regularised boundary matrices M^ e > and N^ e \ This was illustrated 
with a simple example in section 2; in this section we discuss Forman's method in more 
detail and explain how to apply it to the general operator ( |3.10| ). We end the section with 
a summary of the general procedure which we have developed in this paper. 

We suppose, following Forman ||, that the boundary conditions on ( |3.1(J| ) ma y be ex- 
pressed as 





u(a) 




■u(bY 




"0" 


M 


u(a) _ 


+ N 






0_ 



(4.1) 

where M and N are 2r x 2r matrices. This equation is simply the r-dimensional analogue 
of (|2.2|). So for instance, if the boundary conditions are u{a) = u(b) = 0, then 



M 



(4.2) 



I r 0] [0 

, N = 
o oj [I r 

where J r is the r x r identity matrix. 

Now suppose that hi(t); i — 1, . . . , r, is a solution of the homogeneous differential equation 

5Zj Lijhj = 0, and define the 2r x 2r matrix Y(t), which describes the evolution of a solution 

and its first derivative with respect to t, by 

h(a) 

(4.3) 

If y^it), y^it), y 2 (t), are 2r solutions of Lh = 0, (|4.3| ) will apply to each solution sepa- 
rately, that is, H(t) = Y(t)H(a), where 

'y 1 (t) y 2 (t) ... y 2r (ty 

H{t) = (4.4) 

.y x {t) y 2 {t) ••• y 2r (t)_ 

So, in particular, H{b) = Y(b)H(a), or, if the solutions are independent so that detH ^ 0, 



'hit)' 




h(a) 




= Y(t) 




hit). 




h(a) 



Y(b) = H(b)H-\a) 

This explains the second construction (labelled (ii)) in section 2. 

The formula for the ratio of determinants for operators of the type ( |3.10| ) is [|J 

detL e^il^dttrP^P^it)) det(M + NY (b)) 



(4.5) 



(4.6) 

detL exp (| ^dttrP^P^it)) det(M + NY(b)) 

For this result to be applicable, the matrices Pi{t) and Pi(t), and also an d ^2^) 

need not be equal, however the matrix Po{t), multiplying the second derivative, must be 
the same for both operators. In most applications L will be normalised by a L which has a 
different, and simpler, matrix P 2 , but is otherwise the same. In these situations A = Pi, the 
exponential factors in (|4.6| ) cancel out, and the simple formula given by (|2.5|) holds (except, 
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'£(*)' 




~h(b)~ 




= Y(t) 




hit). 




M b ). 



of course, that M,N and Y(b) are now 2r x 2r, not 2x2, matrices). We also note that, 
although the formula seems to be asymmetric with respect to the two points a and b, one 
could just as well define a matrix Y(t) by 

[h(b)~\ 

(4.7) 

so that H(a) = Y(a)H(b). Then det(M+NY(b)) = det(N+MY(a)). Therefore, alternative 
formulae to (2.5) and ( |4.6|) exist, with M and N interchanged and Y(b) replaced by Y(a). 

All of the formalism discussed so far in this section also applies to the problem with 
regularised boundary conditions — the only difference is that M and N are replaced by 
and respectively. We are now in a position to summarise the whole procedure: 

1. Modify the boundary conditions of the original problem by a small amount (e), so that 
yAt) is no longer a zero mode. Let y^'(t) be the eigenfunction of the new problem 
with an eigenvalue \^ which tends to zero as e — > 0. Express the modified boundary 
conditions in the form so that they are characterised by two matrices and 

2. Calculate Y(b) = H{b)H- 1 {a), where H(t) is given by (PQ. 

3. Calculate det(M^ + N&Y(b)). 

4. Calculate from fO). 



5. Hence determine 



lim 



det(M( e ) +N^Y(b)) 
AW 



(4.8) 



6. Calculate the denominator factor det(M + NY(b)) 

7. The ratio of the results of the last two steps gives det' Lj det L. 

We will now study various specific examples where this procedure is applied. 
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V. SPECIFIC EXAMPLES 



The algorithm given at the end of the last section gives a method for determining the 
ratio det'L/ detL. In this section we will give explicit results for a few examples with 
commonly met boundary conditions and also discuss one example in some detail to show 
how the method we have developed works in practice. We will only give results for the 
quantity given by ( |4.8| ), since the final result is found by normalising this by detL, which 
can be found from the formulae given in, for example, Forman's paper M. 

For simplicity we only consider the single component (r = 1) case where the operator 
has the form (|2.1|) , for a variety of boundary conditions. 

(a) With the boundary conditions Au(a) + Bu(a) = 0; Cu{b) + Du{b) = 0, 



det'L 
(Vi\Vi) 



+AC/y 1 (a)y 1 (b), if A,C^ 

-BC/ yi (a) yi (b), ifS,C^0 

-AD/ yi {a) yi (b), if A, £> t^O 

+BD/ yi {a) yi (b), ifB,D^0 



(5.1) 



If all four constants A,B,C,D are non-zero it is easy to see that all four expressions 
are equivalent. Similarly, if only three of the constants are non-zero, then the two 
applicable expressions are equivalent. If only two constants are non-zero, one involved 
in the boundary condition at a, and the other at 6, then only one of the above applies. 
The simple example given in section 2 falls into this class: the boundary conditions 
there correspond to A — 1,B = 0, C = l,D = 0, and in this case ( [BTTl) reduces to 



(b) With periodic boundary conditions u(a) = u{b); u(a) = u{b), 

det'L g/2 (&) ~ 2/2 (a) 
(Vi\vi) yi(a) detiT(a) 

where det H(a) = 2/2( a )yi(o.) — 2/1(0)2/2(0) is the Wronskian. 



(5.2) 
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(c) With anti-periodic boundary conditions u(a) = —u(b); u(a) = —u(b), 

det'L y 2 (b) + y 2 (a) 

(Vi\Vi) Vi i°) detH(a) 



(5.3) 



As an example of the application of these results, we use one of the most well known 
situations where instantons exist: imaginary time quantum mechanics with a potential 
V(x) = \x 2 — ^x 4 Q. As shown in the appendix, this problem leads one to consider 
operators of the form 

d 2 

L = -— + l-3/3 2 dn 2 (dm), (5.4) 
at 2 

where dn is an elliptic function |nj, u = (3(t - t )/y/2 and /? = 2(1 — m)/(2 - m). The 
constants to and m originate from the integration of the second order ordinary differential 
equation which is satisfied by the instanton. The parameter to reflects the breaking of the 
time-translational invariance of the original theory and m is related to the energy of the 
classical particle in the mechanical analogy. The spectral properties of the system can be 



studied by imposing periodic boundary conditions on the path-integral ||11|| , which dictates 
that we use ( |5.2j ) to find the required functional determinant. A straightforward calculation, 
outlined in the appendix, yields 



det'L 2(2 - m) 7 ' 2 
(2/1 1 2/i ) m2 



K(m) E(m) 



2-m 2(1 -m) 



(5.5) 



where K(m) and E(m) are the complete elliptic integrals of the first and second kind re- 
spectively. 

This result simplifies considerably in the limit where the energy of the particle in the 
mechanical energy is zero and consequently the period of the instanton, T, becomes infinite. 
In the appendix it is shown that the asymptotic forms of ( |5.5|) , det(M + NY(b)) and (yi\yi), 
for T large are, respectively, e T /16, — e T and |. Combining all of these results gives 

lim — = (5.6) 

T ^ det L 12 

This is in agreement with previous calculations (e.g. eqn. (29) of Ref. 0). It also illustrates 

the extra complication that may occur if the range (a, b) is infinite. In these cases the 
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numerator (|4.8|) and the denominator det(M + NY(b)) may separately diverge as T = 
(b — a) — > oo. One can avoid these divergences in various ways, but the most obvious way to 
proceed in these cases is to use T as a regulator and to perform all calculations with T large, 
but finite, cancelling out the potential divergences between numerator and denominator 
before taking the T — > oo limit. 



VI. CONCLUSIONS 

In this paper we have developed a simple and effective way of regularising operators 
which have zero modes. The method allows the functional determinants for these kinds 
of operators, with the zero modes extracted, to be calculated. The main advantage of 
the method, and the reason for its power, is that it leaves much of the structure of the 
unregularised problem intact. This means that much of the formalism originally developed 
in this case can be taken over with very little change. The approach which we have adopted 
has not emphasised rigor; it would be very interesting to put this work on a rigorous footing. 
In particular, we have not proved that the results are independent of the precise method of 
regularisation adopted. Until this is done, the results obtained using our method, especially 
for r > 1, should be treated with caution. On the other hand, we have also kept the number 
of examples of the application of the technique to a minimum, preferring instead to give a 
clear and explicit discussion of the methodology. 

Although we have tried to be quite general, describing most aspects of the formalism 
which may arise in practice, there are, inevitably, situations that have not been covered. 
One is the case where there is more than one zero mode present — an example is the model 
studied in 0. The procedure in cases such as this is a simple extension of our previous 
discussion: regularising parameters ei, €2, ■ ■ ■ are introduced for every broken symmetry, and 
hence for every zero mode. This symmetry may be external (spatial or temporal) or internal 
(global or local). The boundary conditions are then modified along the directions of breaking 
by an amount e a , and the prescription given in section IV followed. 
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Our motivation for carrying out this work has been the increasing need to evaluate 
determinants of this kind in many areas of the physical sciences. We hope that the ideas 
presented here are sufficiently straightforward and easily implemented that they will find 
wide application. 
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In this appendix we give details of the calculation of the functional determinant of ( f>. 1| ). 
The motivation for studying an operator of this type is that it arises in the investigation of 
fluctuations about the instanton in one-dimensional quantum mechanics with the potential 



are those of a classical particle of unit mass and energy E moving in the potential —V(x). 
Bounded motion is allowed for E < 0, corresponding to the existence of real instantons. 

Let the values of x at which the particle in this mechanical analogy has zero velocity be 
denoted by a and (3 (0 < a < (3). Then -V(a) = -V(f3) = E which implies a 2 + (3 2 = 2 
and E = —a 2 /3 2 /4. The once integrated equation of motion now reads 
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APPENDIX: 



V(x) = |x 2 — |x 4 ||. The instanton satisfies the equation — 'x + V'(x) = 0, which may be 
integrated once to give \x 2 — V(x) = E, where E is a constant. Solutions to this equation 




(Al) 




— x 



V2 



1 



(A2) 
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where to is the time at which the particle was at x — (3. This may be integrated in terms of 



elliptic functions |T0|: 



x c (t;to,m) = j3dn(u\m) 



(A3) 



where u = (3{t — t )/v^2 and m — 1 — (a 2 / ' (3 2 ). The subscript 'c' denotes "classical" and 
simply indicates that this is a solution of the classical equation of motion 5S/5x(t) = 0, 
where S[x] = J a dt \x 2 + V(x) is the action. The physical significance of the integration 
constant to is clear: since the particle can start at any x (a < x < /3), the time at which 
it reaches j3 (defined to be to) is arbitrary. The constant m, on the other hand, is directly 
related to the energy of the particle, since E = —(1 — m)/2(2 — m) 2 . An alternative to m, 
which also specifies the energy of the particle, is the period T defined by 



T 1 

2 71 



dx 



^J(x 2 -a 2 )((3 2 -x 2 



\/l — m sin 2 9 



m 



1/2 



K(m) 



(A4) 
(A5) 



where K{m) is the complete elliptic integral of the first kind JTC . 

As explained in the main part of the text, we are interested in evaluating the expression 
( p.2| ), and therefore need to determine the values of the functions y\ and 1/2 at the endpoints 
a and b. These two functions are solutions of the homogeneous differential equation Lh = 0, 
where 



L5(t - 1') 



5 2 S 



d 2 

—Tp + 1 - Sx 2 c (t;t ,m) 



S(t - 



(A6) 



Sx(t)5x(t') 

Using the explicit form for x c given by ( |A3| ) we obtain (|5.4|) . But two independent solutions 
of Lh = can be found by differentiating x c with respect to to and m |12|], so we define yi 
and y 2 by: 

dx c (t;t ,m) 



yi(t;t ,m) = 



dt 



(A7) 
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dx c {t;t Q ,m) 

y 2 {t) t , m) = — (A8) 



It is a straightforward exercise in elliptic functions to find from (|A3|) that: 



TflB 2 

yi(t;to,rn) = — y=- sn(w|m) cn(ii|m) (A9) 
yi(t;to,m) = —^—dn(u\m) {cn 2 (w|m) — sn 2 (w|m)j (A10) 

One can now check using (|A5|) that y\{a) = yi(b) and that yi(a) = y\(b) for any initial 
and final times satisfying b — a = T. Therefore, since y% is a solution of Lh = satisfying 
the correct boundary conditions, it is the zero mode for this problem, as expected. 

A slightly longer calculation gives 

dp j dp f3E{u\m) Bu\ 

y2\t; t , m) = - — dn(u m) — < um- -. H > sn( u \m)cn(u \m) 

dm { dm 2(1 — m) 2 J 

f3sv?{u\m) da(u\m) . . 

2(1 - m) ^ ' 

where E{u\m) is the elliptic integral of the second kind. Using the periodicity of the elliptic 

functions 

yi(a) Z \ K[m)m dm 2(l-m) + 2 K[m) S\V2j 

2 (2-m)V 2 \ K(m) E(m) \ 

m \2-m 2(l-m)J 1 ' 

since ft 2 = 2(2 — m)~ l . Here E(m) is the complete elliptic integral of the second kind. 

The Wronskian det H(t) is a constant, and so can be calculated for any convenient t. 

Choosing t = t , which implies u = and so yi(to) = 0, 

det H(t) = y 2 {to) yi(t ) - yi{t ) y 2 (*o) 
f (3 3 m\ ( df3\ 



2 J \dm / 
m 



(A13) 



(2 - m) 3 

Substituting ( |A12| ) and ( |A13| ) into ( p.2[ ), and taking into account the extra minus sign which 



comes about because the operator Q5.4p is minus the definition of operators as given in the 



text, gives ( |5.5|) . 



17 



Following the discussion of the most natural form for L given in section II, we take it 
to be the second functional derivative of the action for the harmonic oscillator with the 
potential V(x) = \x 2 . Then 

d 2 



dt 2 



+ 1 



(A14) 



Choosing yi(t) = e* and ^(0 — e * to be the two independent solutions of the homogeneous 
equation Lh = 0, 

cosh(6 — a) sinh(6 — a) 
sinh(6 — a) cosh(6 — a) 
Using the same periodic boundary conditions which gave (|57 



Y(b) 



(A15) 



det(M + NY(b)) = 2 - Y u (b) - Y 22 {b) 
= 2(1 - cosh T) 
= -4sinh 2 ([2-m] 1/2 K(m)) 



(A16) 
(A17) 



Dividing (|5.5| ) by ( |A17| ) gives the required expression for 

1 det' L 



(A18) 



1 2/i ) detL 

As a check on the results let us look at the limit E — * 0_, i.e. m — > 1 or T — >• oo. In 
this case K(m) ~ | ln(l — m), which from (AD) gives m ~ 1 — 16e~ T . Using E{m) — > 1 as 
m — ► 1, we have 



det^L 
(2/i 1 2/i ) 



— , as T -> 00 
lb 



(A19) 



Since from ( |A16| ), detL 



-e as T — > 00, 



lim 



1 det' L 



1 
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(A20) 



T ^°° (2/1 1 2/i) detL 

The sign is the expected one: the zero mode which has been extracted, y\, has a single 
node, which leads us to deduce that L has only one eigenfunction with a negative eigenvalue; 
all the other eigenvalues are non-negative. The signs of (|A17 ) and ( |A19|) are not those that 



we might naively expect, but these signs have no meaning separately — both the magnitude 
and sign of these terms can be changed at will by the replacement M — > AM, N — > XN, 
where A is any real number. 

The ratio ( |A20| ) agrees with the calculation of Ref. |J. To see this we note that a — > 
0, /3 — > v^2 as m — > 1, hence the instanton becomes 



x c (t; t Q , m = 1) = \pl sech(t — t ) 

yi(t; t , m = 1) = v2 sech(t — t ) tanh(t — t 



lim ( yi \ yi ) = - 

1 — >oc J 



(A21) 
(A22) 
(A23) 



Combining (|A20|) and (|A23| ) gives ( p.6| ), as required 
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